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Abstract 



In the standard picture, the inflationary universe is in a supercooled state 
which ends with a short time, large scale reheating period, after which the 
universe goes into a radiation dominated stage. An alternative is proposed 
here in which the radiation energy density smoothly decreases all during an 
\ inflation-like stage and with no discontinuity enters the subsequent radiation 

dominated stage. The scale factor is calculated from standard Friedmann 
■ cosmology in the presence of both radiation and vacuum energy density. A 

large class of solutions confirm the above identified regime of non-reheating 
inflation-like behavior for observationally consistent expansion factors and 
not too large a drop in the radiation energy density. One dynamical realiza- 
Q^' tion of such inflation without reheating is from warm inflation type scenarios. 

O ■ However the solutions found here are properties of the Einstein equations 

with generality beyond slow-roll inflation scenarios. The solutions also can 
be continuously interpolated from the non-reheating type behavior to the 
standard supercooled limit of exponential expansion, thus giving all interme- 
diate inflation-like behavior between these two extremes. The temperature of 
the universe and the expansion factor are calculated for various cases. Im- 
plications for baryongenesis are discussed. This non-reheating, inflation-like 
regime also appears to have some natural features for a universe that is be- 
tween nearly flat and open. 
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I. INTRODUCTION 

In the original conception of inflation |IJ, it was assumed that the universe underwent 
isentropic expansion during the stage of rapid growth of the scale factor. The entropy 
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required to make the post-inflationary universe consistent with observation was assumed 
to be generated in a short-time reheating period. However, it is clear that for a range of 
moderate thermodynamic conditions, the cosmological horizon and flatness problems, which 
are explained by inflation, require only the kinematic property that the scale factor grows 
rapidly. More recently it was realized that these kinematic conditions could still arise 
in the presence of a sustained radiation component during inflation. Specifically, in [0] it 
was shown that under certain isothermal conditions inflation could still occur. More so it 
was shown there, that within these limits, the initial seeds of density perturbations could 
be dominantly of thermal instead of quantum origin. A realization of an isothermal or 
warm inflation scenario, in the context of slow-roll scalar field dynamics for parametrically 
large dissipation, was shown in || to be consistent with observational constraints for the 
amplitude and expansion factor, without requiring a ultra-flat Coleman- Weinberg potential, 
which in order to form requires the coupling constant to be fine tuned. Questions about the 
fundamental origin of large dissipation are still left open. 

The warm inflation scenario served as a explicit demonstration of an otherwise true but 
ambiguous statement, that inflation can occur in the presence of a thermal component. That 
this is true is self-evident, as for example within the context of scalar field theory. Here the 
requirements for exponential expansion are 

pv > Spt (1) 

and 



Pv 3> Pkinetic (2) 

with 

Pkinetic = + Pv (3) 

Here bp^ is the energy density perturbation, and p v and p r are the background vacuum and 
radiation energy densities respectively. Thus energetics alone does not prohibit the relation 

Pv > Pr > ^0 2 , Sptf,. (4) 

By itself this inequality gives no indication on the extent that radiation can modify the 
supercooled scenario. However, the warm inflation scenario || demonstrated that at least 
in the limit of near thermal equilibrium, the effect is nontrivial. By reexamining this scenario 
solely in terms of energetics in |4]], it became evident that both supercooled and thermal 
slow-roll scenarios could be viewed as limiting cases of a class of nonequilibrium kinetic 
possibilities. A preliminary step to a nonequilibrium study is determining the possible 
kinematic behaviors of the scale factor for a universe in a mixed state of radiation and 
vacuum energy. This is the first motivation that leads us to examine the scale factor in this 
paper. 

In light of this, we find it useful to distinguish between the behavior of the scale factor, 
which we consider kinematics, from the underlying dynamics that induces this behavior. 
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The classification of scale factor behavior is considered kinematic, because it involves char- 
acterizing different solutions and different regimes of a given solution, all arising from a 
particular equation. Besides inflation, common amongst these are radiation dominated and 
matter dominated behavior. Originally inflation was associated with an exponentially grow- 
ing scale factor Subsequently any form of accelerated expansion (R(t) > 0) has become 
associated with inflation. 

Dynamics enters in determining the time evolution of the background stress energy ten- 
sor, which is the driving source in the scale factor equation. In the context of dynamics, 
inflation may or may not arise due to a phase transition. In general, dynamics is stochastic, 
although the degree of stochasticity may well be approximated by pure dynamics or near- 
equilibrium statistical dynamics Q. Inflation scenarios realized in a supercooled regime are 
examples of the former, whereas warm inflation scenarios || are examples of the latter. 

The present most successful formulation of supercooled scenarios is new inflation |5]||. 
Although several variants of the original scenario have been formulated (for a review please 
see |7|||), U P to observation the basic assumptions and mechanism are the same. The 
new inflation assumptions are that dynamics can be described by a suitable potential with 
a suitable order parameter, known as the inflaton, and that evolution is governed by the 
Lagrangian equations of motion. The basic mechanism of new inflation is slow-roll dynamics 
at supercooled temperatures. 

In the simplest form of new inflation, the inflaton is a scalar field. The conventional 
treatment of scalar field dynamics assumes that it is pure vacuum energy dominated. The 
various kinematic outcomes are a result of specially chosen Lagrangians. In most cases the 
Lagrangian is unmotivated from particle phenomenology. Clear exceptions are the Coleman- 
Weinberg potential with an untuned coupling constant, which is motivated by grand unified 
theories ||, and supersymmetric potentials, although in the latter case, the choice of the 
supersymmetric potential is again arbitrary, and in the former case new inflation is incon- 
sistent with observation. Making one extension to the new inflation picture, the behavior of 
the scale factor can also be altered for any given potential when radiation energy is present. 
Out of pure kinematic interest, this effect has reason to be examined. 

More so than just this reason, one may also project to circumstances sometime in the 
future when observational data will allow determination of the optimal potential amongst 
the candidate choices (for examples of recent attempts please see |jn,0]). If one accepts the 
new inflation approximation that the relaxational dynamics of the inflaton can be described 
by a potential, the next question is what is the microscopic origin of this so preferred 
potential. If one were restricting to supercooled scenarios, one argument is that the so 
preferred potential happens to be the one that formed during the rapid quench at the onset 
of inflation. Another argument is that this is a fundamental zero temperature potential 
of an elementary field in the Lagrangian. Since for supercooled new inflation scenarios, 
one of the unanswered questions is that no potential that is suitable for inflation has an 
already known phenomenological origin, the second argument is highly predictive. Yet to 



x The role of stochasticity in cosmology has been emphasized by the Maryland school. For a review 
please see pUp^ ] 
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substantiate either claim, one would need to study the evolution of the potential from its 
high temperature state during the quench. This would lead to examining the interplay 
between radiation and vacuum energy density at the onset of inflation. 

Having appreciated this point, the time interval in which this transition occurs becomes 
important. The short time regime is relevant to supercooled new inflation scenarios and the 
extent to which this interval can be extended is relevant to warm inflation scenarios. Thus, 
whether stated in the conventional sense of new inflation or the extended sense of warm 
inflation, the out-of-equilibrium evolution of the inflationary potential will require study, 
and as an initial step, the scale factor dynamics needs to be examined in a mixed state of 
background vacuum and radiation energy density. This additional connection to supercooled 
scenarios provides a second motivation for this study. 

To completely study the nonequilibrium dynamics, the problem divides into two steps. 
The first step is determining the regimes in which accelerated expansion and pure inflation 
can occur and characterizing the behavior of the scale factor in these regimes. The second 
step is understanding within the allowed regimes, the class of spectra of primeval energy 
density perturbations. The first step is moderately model dependent and mainly involves 
energetics and Friedman cosmology. The second step is a more acute problem of dynamics. 
Although we will only address the first step in this paper, let us make a few comments about 
the second step. 

In general there is no unique formulation of nonequilibrium dynamics for almost any 
system. The first step in formulating any approach requires understanding the scales in 
one's problem. For inflation the simplest assumption is that there are two scales: a longtime, 
long-distance scale associated with vacuum energy dynamics and a single short-time, short- 
distance scale associated with a random force component. The Hubble time during inflation, 
l/H, appropriately separates the two regimes. For grand unified theory || , this time interval 
is 1/H as I(T 34 sec. 

The assumption of a longtime scale for the evolution of the vacuum is based on obser- 
vation. Otherwise inflation would not sustain itself sufficiently long nor would the energy 
release maintain smoothness. Accepting this as an empirical constraint, the relaxational 
dynamics of the inflaton's order parameter justifiably could be described by a free-energy 
functional. What the specific functional is requires dynamics. In the presence of a radiation 
component, the functional need not have any similarity to a fundamental potential from the 
underlying quantum field theory. Furthermore, in grand unified theories as an example, the 
characteristic time scale of inflation, 1/H, is about 10 10 times faster than the characteris- 
tic hadronic interaction scale (1/Aqcd), which is a comparison scale where there is good 
empirical understanding about matter. Thus at the inflation scale, familiar concepts about 
matter and from field theory about near-thermal-equilibrium-motivated effective potentials 
also need not be appropriate. 

The problem here has similarities to certain phase separation problems commonly known 
in association with binary alloys, and the name synonymous with them, spinodal decompo- 
sition, has been used before in new inflation cosmology PL The similarity in both cases is 



2 For scalar field inflaton dynamics, the analogy actually is to spinodal decomposition for a non- 
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that the system is being cooled faster than its characteristic response time to equilibrate. 
Of course if this analogy is meant to be complete, the cooled system should still be at a 
nonnegligible temperature, since at least in the binary alloy problem, the relaxational dy- 
namics is driven by short ranged thermally excited fluctuations. The analogy to spinodal 
decomposition not only gives a nice guiding picture, but it also has a type of consoling 
appeal, which covers for our ignorance about matter, much less quantum field theory, under 
such extreme conditions, since at least in the context of alloys, the problem is considered 
sufficiently complex to make phenomenological modeling of the nonequilibrium potential an 
accepted practice. If viewed in the same way, the several scalar inflationary potentials that 
have been suggested could be interpreted as the cosmologist's attempt at nonequilibrium 
phenomenology. 

A. Hypothesis 

Although the reasons given above well motivate examination of the scale factor, I will now 
describe an alternative to the standard inflationary universe scenario. Consider the following 
possibility which will be demonstrated in the sequel. It should be easy to convince oneself 
that a radiation energy density p r (t) of say one part in ten thousand to the vacuum energy 
density p v (t), probably should not alter too much the inflation- like behavior of the scale 
factor. However when looked upon in terms of the temperature T r of the radiation energy 
density, this implies that T r is only an order of magnitude below the scale of the vacuum 
energy density. If such a state for the radiation energy density could be maintained by the 
mutual effects of constant vacuum energy decay and a steadily decreasing acceleration of the 
scale factor, it could be possible for an inflation- like stage to smoothly enter into a radiation 
dominated stage without any discontinuities in p r (t). This possibility was suggestive from 
formulating the warm inflation scenario [[3]. 

In this paper evidence is presented for inflation-like trajectories of the scale factor which 
solve the horizon and flatness problems, but for which the radiation energy density monoton- 
ically enters the post-inflation radiation dominated stage with in particular no intermediate 
reheating stage. This is a regime in between the radiation dominated and inflation regimes, 
which has features similar to both a big-bang-like explosion and an inflation-like expansion. 

The paper is organized as follows. In the next section the problem is formulated, general 
solutions are given in section 3, special examples are given in section 4 and finally the conclu- 
sion is in section 5. From the class of solutions that we find for the scale factor, supercooled 
expansion, which we call supercooled inflation, is a limiting case. This is a kinematic iden- 
tification. A particular and most noteworthy dynamic realization of supercooled inflation 
scenarios is the class of new inflation scenarios. Supercooled inflation has associated with 
it also a range of power law [12], quasi-exponential [13] and exponential |5||| behavior for 



the scale factor. However these varied behavior arise from the specifics of the particular 



conserved order parameter, such as found in certain domain growth problems 1 33 , 34 ] , whereas the 
binary alloy problem involves a conserved order parameter 
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Lagrangian that is being considered. In what we will examine, for any given Lagrangian, a 
large range of behavior may still arise, depending on the radiation energy density content. 

There is one other problem that the present work may help clarify. We will discuss it 
briefly here. However it gets into the realm of field theory dynamics, which this paper will 
mainly avoid. The most notable shortcoming of new inflation is in explaining small scale 
energy density inhomogeneities |T^j6|JT5HT7| . The problem is sometimes referred to as the 
amplitude fluctuation problem |I5|. The warm inflation scenario in |3j is a solution to this 
problem. However our formulation there did not detail a time history for an inflation-like 
state with radiation. The present work does and in fact was its starting motivation. However 
due to the generality of the solutions given here, it appears better to consider warm inflation 
as a particular dynamic realization within the big-bang-like inflation regime. 

Our equations can also be examined for the initial stage of entering into the rapid ex- 
pansion state, but we will not study that here. 



II. FORMULATION 



We are interested in the scale factor from some short time after the initial singularity, 
when quantum gravitational effects become negligible. We assume space is homogeneous and 
isotropic, and restrict ourselves to Friedman cosmology with the Robert son- Walker metric 



ds 2 = dt 2 -R 2 (t) 



dr 2 



1 — At 2 



+ r 2 d6 2 + r 2 sin 2 



(5) 



For notational convenience, the origin of cosmic time is defined as the beginning of our 
treatment. 

Let us start with the standard equations of Friedman cosmology for the scale 

factor R(t) in the presence of vacuum energy density p v (t) and radiation energy density 
p r (t). The equations of state which relate the energy density p to the pressure p are 



p v (t) = -p v (t) 



(6) 



Pr(t) = -Pr(t). 



(?) 



In Friedman cosmology the ten Einstein equations G^ v = RnGT^ reduce to two indepen- 
dent ones, which are from the time-time component, also known as Friedmann's equation 



R 2 k 8nG 



and from any of the three diagonal space-space components, all of which give 

R R 2 k 
2 R + V + 1V=-^ GP - 



(8) 



(9) 



For our purposes it is preferable to use two other equations obtained from these, the scale 
factor equation 
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R 8nG 



R 3 

and the stress energy conservation equation 



\p v {t) - p r {t)\ , 



R(t) 

p r (t) = -4p r {t)-^--p v (t), 



(10) 



(11) 



where we have used the equations of state eqs. (|) and (|7|). We aim to solve for R(t) and 
p r {t) in eqs. (pi]) and (|TT1), for a prescribed p v (t), for t > 0, and with arbitrary initial 
conditions for R(t), R(t), and p r {t) up to the constraints 



and 



R(0) > 



R(0) > 



p r (0) > 0. 



(12) 
(13) 



(14) 



By taking the sum and difference of Friedman's equation eq. (§) and the scale factor 
equation eq. (|T0[) , the vacuum and radiation energy densities can be separately expressed 



in terms of the scale factor as [TB5 



Pv(t) 



16nG 



R & Jt_ 
R + R2 + R2 



(15) 



Pr{t) 



16nG 



_R R 2 k 
~R + R 2 ~ + R 2 



(16) 



For an arbitrary test vacuum function p v (t), one can use eq. flT5"D to solve for R(t). 
We make the substitution 

s(t) = R\t). (17) 

Eq. (p~5|) then becomes the inhomogeneous wave equation with time dependent frequency 

32ttG 



s — 



-p v (t)s = —2k. 



(18) 



This equation has been widely studied JSyflll. Again we are interested in the solutions to 
eq. (18) for t > with arbitrary initial conditions for s(t) up to the constraints from eqs. 
(0)-(rrj) and (|17D which imply 



and 



s(0) > 
s(0) > 0. 



(19) 
(20) 
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III. SOLUTION 



In this section solutions are obtained for the scale factor from eqs. (|i~0D and (|TT| ) for a 
large class of vacuum energy decay functions. Even before getting this specific, there are 
two general features, one at short and one at long time, which are recurring themes to the 
existence of the big-bang-like inflation regime. At long time, if p v (t) goes to zero sufficiently 
fast, from eq. (jTSf) and eq. ([17]) one can see that R(t — ► oo) ~ t 1 / 2 , thus tending to a 
radiation dominated behavior. At short time, for an initially radiation dominated universe, 

Pr {t ~ 0) > Pv {t ~ 0), (21) 

eqs. flTOD and (|T]]) imply R(t ~ 0) ~ (a + bt) 1 ^ 2 . Alternatively, this can be seen from 
Meissner's separation eqs. (]T5|) and (0), since eq. (]2lf) implies from eqs. fllSf) and (p~6f) that 

s 2 (t ~ 0) > s(t ~ 0). (22) 



Taylor expanding s(t) about the origin as s = So + sit + S2t /2 H — • s n t n /n\ H — • +, eq. (|22|) 
implies s 2 3> S2. Using this and eq. (pT8|), one can study the initial condition dependence of 
entering the inflation- like stage, but we will not pursue that here. 

Let us now turn to specific solutions. Although this paper is focused on the kinematic 
possibilities for the scale factor, independent of justification from any specific field theory, we 
will motivate a class of vacuum decay functions from a general class of scalar field dynamics. 
In fact as we will show below, in the limit of strong dissipation, this motivation can be partly 
justified. 

We consider stochastic evolution for the inflaton governed by the Langevin-like equation 



<j>(t)+V(<f>(t))=r)(t), (23) 



R(t)_ 

where rj(t) is a random force function with vanishing ensemble averaged expectation value 

(rj(t)) = 0. (24) 

The effect of the inflaton's interaction with radiation is represented by the dissipative con- 
stant T and the random force function t](t). In a simple model for the radiation system and 
in the limit of pure inflation (R/R =const.) this equation was obtained from quantum field 
theory in M. 

We are interested in the limit of strong dissipation 

r » | (25) 

and the slow-roll regime 

m » h- (26) 

For our present purposes, the ensemble averaged equation of motion is all that we need. 
Thus in the above specified limits, eq. (1231) becomes 
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~dt 



1 dV(4>) 



Let us consider potentials of the form 

V(<f>) = AM 4_n (M 

in the region 

< 6< M 



(27) 



(28) 



(29) 



where A is dimensionless. For inflation driven dynamically at the grand unified scale M ~ 
M GUT « 10 14 GeV. 

Globally all of the potentials in eq. (p8j ) are improper for slow-roll inflation scenarios, 
since they fail to represent symmetry breaking. However our present interest is the behavior 
of the scale factor for a large class of slow-roll conditions. In this sense eq. ( p8|) represents 
a class of local approximating potentials from which an arbitrary potential can be piecewise 
constructed. Thus it is also not a concern that such potentials have no minima for odd n 
and are nonanalytic for noninteger n when <fr = M. 

In fact near the global minima, where the vacuum energy goes to zero, quadratic (n = 2) 
dependence would be the normal expectation for any generic free energy functional. This 
case is not only of special physical interest but is also mathematically a little different. We 
will differentiate this case of n = 2 from all others and refer to it as the quadratic limit. 

To keep our discussion explicit, we will express the results that follow in the context of 
the slow-roll inflation scenario. However it should be noted that the solutions for the scale 
factor given below carry a relevance beyond the slow-roll scenario. Let us briefly recall the 
slow-roll scenario. In the standard setting of the slow-roll transition, the inflaton starts near 
the origin and is making its decent to the symmetry broken minima at <fi = M. 

At the origin of cosmic time we will assume the slow-roll transition begins with 

0(0) = eM (30) 

and e<l. With these initial conditions, the solutions of eq. ( |2~TD for potentials in eq. ( |28|) 

are for n = 2 



<j){t) = M 



and for n ^ 2 



with 



4>(t)=M l-(B n (t + t 0n )) — 



B, = 



4AM 2 



(31) 



(32) 



(33) 



and 



9 



B. n = 



n(n - 2)AM 2 



(34) 



Here t 02 and to n are suitably adjusted to satisfy eq. fl3"0"|) . Equating the potential to the 
vacuum energy density 



implies for n = 2 



and for n ^ 2 



p v {t) = V(<f>(t)) 



p v {t) = XM 4 exp(-B 2 (t-t 02 )) 



p v (t)=\M 4 [B n (t + t 0n )}^ . 



(35) 



(36) 



(37) 



Substituting the above in eq. ([18]) and solving the homogeneous (flat space) equation, we 
find for the scale factor from eq. (|17D for n = 2 



R = ^C 1 I (z 2 (t)) + C 2 K (z 2 (t)) 
and in the n^2 case for all but n = 4 



where 



R= [B n {t + t 0n ))* JCJ^iznit)) + C 2 K^(z n (t)), 



(38) 



(39) 



(40) 



with 



= 4(2 - n)H n t 0n t_ \ — 
n{ ) ~ (4 - n) W + V 



8ttGAM 4 (B 2 t 02 
H 2 = \ exp 



(41) 



(42) 



/ 87rGXM 4 

R n = \ o (5„t J 2(2 -"'. 



(43) 



In eqs. 
definition 



and ([43]) we have identified the Hubble parameter at t = based on the 



H 



(44) 
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for the respective vacuum energy densities in eqs. (|36|) and fl37|). In the Appendix we have 
listed properties of Modified Bessel functions that will be useful to us. Finally from the 
cases for n = 4 the solution is 



R = \B A {t + to^yC^t + toJl* + C 2 \B 4 {t + * 04 )|-" (45) 

where 




128vrGAM 4X 1/2 



A*=o(l+ \m ■ ( 46 ) 



The inhomogeneous wave equation in eq. (|TB|), which is for curved space k ^ 0, can be 



solved from the above solutions for the homogeneous equation by familiar methods pC||2T 
Irrespective of the slow-roll scenario, the results eqs. (^) and (|39"D are valid for any scenario 
that motivates vacuum decay behavior as in eqs. (^) and fl37|). Likewise for other types of 
vacuum decay functions, eq. ( |i~8D can be solved |. 



In the next section, the solutions eqs. (|38|) and fl39|) will be studied through specific 



examples. Here some of their general features will be noted. The quadratic limit is examined 
first. The growing mode as t — > oo in eq. (|38|) from eqs. ( |A1| ) and (|A2| ) is K (z 2 (t)) with 



R{t -> oo) ~ £ 1/2 , (47) 

thus asymptotically exhibiting radiation dominated behavior. 

Inflation-like expansion at intermediate time is also governed by K (z 2 (t)). From eq. 



( A~2|) a large expansion factor of e N with AT > 50 will require 

2H 2 



Bo 



N. (48) 



This is like what one would expect, since the vacuum energy density must decay sufficiently 
slowly relative to the expansion time for the scale factor, in order to be the driving source 
for inflation- like behavior in the Einstein equations. 

In order to establish the dominance of the ^0(^2(0) term to the Io{z 2 (t)) term in eq. 
(|58D, what remains is to show that there is no way for the initial conditions to force C\ 
to be exponentially large relative to C 2 - Treating 4H 2 /B 2 ^> 1 and using eqs. (|A5|) and 
this follows from the constraints eqs. ( |i~9|) and (|20|). Eq. (^) could be satisfied for 



C± exponentially large but negative relative to C 2 , but then eq. (|19D would not be satisfied. 



One case is during reheating in supercooled scenarios. For this, the vacuum decay function in eq. 
(|j~8| ) should have the approximate time dependence e~' r '*(l + cos Bt) with B ~ M S> H. These 
types of equations are treated in pC| . The |T| = case is the Mathieu equation. The solutions 
of these equations describe the scale factor behavior during those stages of reheating when the 
equation of state eq. ([]) is valid for the inflaton. Another nontrivial aspect of scale factor behavior 
in supercooled scenarios is at the beginning where initial condition dependence on pre-inflationary 
radiation energy density can be studied. For this, the solutions eq. ( p9| ) for n ~ are useful. 
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Having established the dominance of the Ko(z2(t)) mode, let us estimate the expansion 
factor for a single n = 2 section of potential eq. (|2lf) with e = in eq. ( |30"1) so that t 02 = 
in eq. fl3~T|). We find for the asymptotic behavior 

R(t - oo) (£ 2 t/2) 1/2 e 2W 

i?(0) ~ (7rB 2 /8H 2 y/* 1 J 

so an expansion factor of e 2H2 ^ B2 . 

Away from the quadratic limit (n 7^ 2) from eqs. (|39|), (|A5|) , (|A6|) we see that for 

^—^ > (50) 
2-7i v ; 

the solution will grow exponentially at large time, thus never asymptotes into a radiation 
dominated behavior. This corresponds to a vacuum decay function that decays slower than 
1/t 2 at large times in eq. fll8l) . Radiation dominated behavior at large time is attained for 

4 -"<0 (51) 



2-n 

which implies potentials in eq. (|28[) with 

2 < n < 4 (52) 

or vacuum decay functions in eq. ([0|) that decay faster than 1/t 2 . Finally for n = 4, which 
corresponds to a vacuum decay falling-off exactly as 1/t 2 , R(t) has the same power law 
behavior throughout, with a growth bounded from below by t 1//2 . As such, this case is not 
useful for our present purpose. This implies that the only symmetric potential about the 
symmetry broken point, <fi = M, that leads to radiation dominated and not inflation-like 
asymptotic behavior is the quadratic case n = 2. As an aside, note that the n = 4 case 
is interesting since on either side are solutions with two very different types of asymptotic 
behavior 0. 

Returning to the cases in eq. fl52|) , the growing mode in eq. (p9|) is if 2-n (z n (t)). Let us 

— — 4 ~ n 

estimate the expansion factor for a single n^2 sector of the potential eq. fl28|) in the range 

eq. (|52l) for e = in eq. ([$(]) so that 



B n t 0n = 1 (53) 
in eq. Cj32|) . The arguments are the same as above for the n = 2 case with the final result 



4 It is also an interesting coincidence that n = 4 separates renormalizable and nonrenormalizable 
scalar quantum field theory, with the nonrenormalizable side, n > 4, corresponding to the observa- 
tionally inconsistent non-radiation dominated asymptotic scale factor behavior. Furthermore the 
n = 4 case neither asymptotes to radiation dominated behavior nor is believed to be nonpertur- 
batively a nontrivial quantum field theory [35]. Of course for the inflaton, since it is coupled to 
gravity, the whole theory is always nonrenormalizable in any case. 
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R(t 



oo 



R(0) 



{B n t) 1/2 exp 



'2{n - 2)H n 
. (4 - n)B n 



(54) 



Before closing this section, one additional qualifying statement is needed about the solu- 
tions eqs. ( |3g|) and fl39[) if the vacuum decay functions eqs (|36|) and (^) are obtained from 
slow-roll scalar field dynamics. Recall that the energy density and pressure of the zero mode 
of the scalar field are 



1 



P4> = of + V(0) 



(55) 



P0 = ^0 2 -^(0). (56) 

Therefore, the equation of state eq. (|]) is valid in the limit that the potential energy 
dominates the kinetic energy 

\^ « Vffl. (57) 

One must check that this kinetic energy suppression condition is always valid. 

For this, first recall that the exact equation of motion for the inflaton in the limit eq. 
(p5p is the second order equation 

$ + T0 + V'(4>) = 0. (58) 

For the quadratic case n = 2, this equation can be exactly solved and it can be verified that 
the slow-roll condition eq. fl2"6j ) and the kinetic energy suppression condition eq. ( |5"7D are 
both valid for all t > provided 

AM 2 , s 

-jT" « °- ( 59 ) 



In addition, if p v (0) is required to be large in eqs. (|36|) and (|37|), so that A can not be made 
tiny, eq. ( p9| ) implies 

r > M. (60) 

This is the large dissipative regime required for warm inflation 0]. 

For the cases 2 < n < 4, to verify eq. fl5"7|), first it will be shown that the solutions eq. ( |3"2"D 
are consistent with the slow-roll condition eq. ( p6|) for all i > 0. Next a direct verification 
will be made that the solutions eq. ( |3~2"1) respect the condition eq. (|57|) . Addressing step one, 
it is observed from eq. (|32"D that for the entire range 2 < n < 4, </>(t) vanishes faster than 



4>{t) as t — » oo. Thus eq. (^3) is satisfied under the same parametric restrictions as in the 
n = 2 case, eqs. (|5Tf ) and (|6ll|). Proceeding to the second step, it can be verified from the 
slow-roll approximate solutions eq. ( j3~2"l) that 4> 2 {t) vanishes faster than V(<fi(t)) as t — ► 00. 
Thus in the regime eq. (j60j) , eq. (|5T| ) is satisfied for all t > so that eq. (|) is always valid. 

To summarize, it has been verified that the equation of state eq. (|6|) is valid for the 
scalar field for all t > and in the entire range 2 < n < 4, when in the strong dissipative 
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regime eq. (|60|). This type of slow-roll motion is analogous to an over-damped oscillator [|J. 
Note that confirming the validity of eq. (§) for all t > is more than needed, since in any 
case p r (t) overtakes p v (t) at some much earlier stage. 

The results presented in this section now demonstrate the existence of inflation-like 
scale factor trajectories which smoothly go into a radiation dominated behavior without a 
discontinuous reheating stage. 

IV. EXAMPLES 

In this section we will examine some specific examples from the solutions for the scale 
factor in eqs. (|38l) and (^). In these examples we will see how the radiation energy density 
eventually overtakes the vacuum energy density with no discontinuities, and in the processes 
the universe smoothly goes from an inflation- like to a radiation dominated stage. We will 
also study the magnitude of decrease in the radiation energy density, thus the temperature 
of the universe, from before to after the inflation-like stage. In supercooled scenarios, the 
post-inflation temperature is referred to as the reheating temperature, but here it is better 
to call it the initial temperature after inflation, Tai- In particular the inflation-like stage is 
defined as the time period when the scale factor has positive acceleration 

R{t) > 0, (61) 

with the time "just before" , t BI , and "just after" , t^i, inflation being defined as the endpoints 
of the accelerated expansion interval 

In subsection 4b we will examine a particular n ^ 2 case from eq. ([39]) which can be 
fully expressed with simple analytic functions. Then in subsection 4c we will examine the 
quadratic limit. For this study, we will first convert to a set of dimensionless quantities. 

A. The Dimensionless Theory 
We will work with the dimensionless quantities defined as 

a(T) s R^T) (62) 

c(r) . #^ (64) 



where dimensionless time 



Pv{tbi) 



t = Ht, (65) 
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tbi is the time when accelerated expansion begins, and H is defined in eq. (f44|) except with 
the vacuum energy density evaluated at t bi , p v {j BI ). Defining s a {r) = a 2 (r), the Meissner 
separation, eqs. ( |T5"D and (]T6|), in terms of s a (t) and the rescaled quantities is 

, / \ 1 d 2 s a {r) k , . 



and 



2 



/ \ _ 1 d 2 s a (r) 1 / ds a (r) . , > 

lJ 4 Sa (r) rfr 2 4sI(t) { dr ) + 2HH a {rY 1 J 

The radiation energy density will be related to a temperature measure by the Stefan- 
Boltzmann radiation law 

Pt (t)~T\t). (68) 

This law need not hold under far from equilibrium conditions, but we will nevertheless refer 
to T(r) as the temperature of the universe at time r. We will study the temperature of the 
universe in terms of the ratio 

a(r) = (69) 

The field theory quantities will also be rescaled. The natural scale for them is M not H 
and in general these two scales are different. This scale disparity is an inherent feature of 
scalar field slow-roll dynamics. A primitive source of the dilemmas encountered in slow-roll 
scenarios is its two scale nature. The natural time scale in the field theory to release the 
vacuum energy 1/M, in general differs from the characteristic cosmological expansion time 
1/H. In grand unified theories this disparity works against theoretical preference, since 

1/M < 1/H. (70) 

Had this inequality been reversed, it would have been parametrically satisfying and per- 
haps a strong argument for theoretical consistency between cosmology and particle physics. 
However since this is not the case, it either means slow-roll dynamics is wrong, field theory 
dynamics for inflation at the grand unified scale is wrong, grand unified theory is incomplete 
or wrong or that the physics needs further elaboration, perhaps from nonequilibrium meth- 
ods. We will not address the dynamical problem here, but it is worthwhile to keep track 
of the scale disparity. Thus we will rescale everything with respect to H, but for quantities 
where M is the natural scale, the rescaling will include the additional factor 

The field theory quantities are rescaled as 

T = 1 (5H (72) 

and 

(j) = af3H. (73) 
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B. n=8/3 



Let us consider the case n = 8/3 from the n ^ 2 class of potentials in eq. (2£j) which in 
rescaled parameters is 



V{a) = AM 4 (1 - <x) 8/3 . 
Solving the slow-roll equation of motion 



where using eqs. (|3 



«8/3 



and with the initial condition 



we find 



_ ^8/3 = 16A/3 

Fg/3 97 
(7(0) = 0, 

1 

( K8/3 r + 1)3/2- 
This implies that the rescaled vacuum energy density is 

Kr) = 7 T7- 

1 ' (k 8/3 t + 1) 4 

Solving the homogeneous wave equation eq. ( |I8"D using eq. ([79] ) we find 



a(r) = 1 



or 



\ 



(k 8/3 t + 1) 



Ci exp 



2r 



C 2 exp 



-2r 



(« 8/3 r + 1) 



(k 8/3 t + 1) 

Here tbi, Ci and C2 are determined by the initial radiation energy density r, 

c(0) _ 



the defining relation for re/, 



and from eqs. (|62| ) which implies 



6(0) 



6(rs/) = c(tb/) 



a(r B i) = 1. 



(74) 



(75) 



(76) 



(77) 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



From eqs. (^g) and (|67| ) and restricting to flat space, explicitly the first two conditions above 
imply from eq. fl31~|) 
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'ds a (o)\ 2 



and from eq. 



4,2(0) V dr 



'ds a (r B i, 



r + 1 (84) 



% s K t bi) V dr 



2 

I =b(r BI ). (85) 



In eq. ( jS5|) s 2 a {jBi) = 1, but we retain it explicitly since the same equality holds at tai, and 
we will use it below to determine tai- 

Let us verify the various general features discussed in earlier sections for this specific 
example. At large time, by inspection of eq. (|80| ) one finds 

a(r -> oo) ~ r 1/2 , (86) 

which verifies an asymptotic radiation dominated behavior. Turning to the growth of the 
scale factor, to obtain an exponentially large one from eq. (|SU[) at large time, the only way 
is from the first term (the growing mode) and only if l//Cg/3 ^ 1- The constraints eqs. 
(|T9|) and fl20p imply C\ > |C 2 |, so that the growing mode can not be suppressed due to an 
exponentially small coefficient C\ relative C<i. To leave no ambiguity, let us focus on one 
case amongst a large class which are all about the same for what we want to study. For 
simplicity, in the case we will consider, arrange the initial conditions so that inflation begins 
at the origin of cosmic time tbi = 0. One can confirm from eqs. (|83|) and fl34|) that C\ and 
C*2 are about the same order of magnitude, so that at large time, only the growing mode 
need be retained. For completeness we find 

Cl = \(l + V2- K ^) (87) 

and 

ifl-V2+^V (88) 



C 2 

2 

From the growing mode in eq. (|80|) we see that to obtain N e e-folds of expansion requires 

— = N e , (89) 

^8/3 

which agrees with our general n ^ 2 approximation formula eq. ([54]). In grand unified 
theories, for the Coleman- Weinberg potential with an untuned coupling constant, one has 
(3 ~ 10 4 and A ~ 1 so that from eqs. (|76|) and (|59|) this implies 

7 ~ A^ e 10 4 , (90) 

which for N e ~ 50 implies 7 ~ 10 6 . This gives an overviewed explanation for the large 
dissipative constant found in the warm inflation scenario of 0. However the estimate here 
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for 7 is a little higher, because in the actual scenario, the finite temperature Coleman- 
Weinberg potential has a smaller curvature. In this simplified discussion, this means (3 is 
smaller. 

We see once again that the largeness of (3 makes a seemingly undesired appearance in the 
dynamics. Whereas in supercooled scenarios it forces a fine tuning of the coupling constant, 
here it forces the dissipative constant to be larger than one would naively want. It remains a 
theoretical question whether such a large dissipative constant has an explanation. However 
up to the parametric level, the warm inflation approach still has a hope of salvaging the 
simplest grand unified theory motivated Coleman- Weinberg model for inflaton dynamics. 
In new inflation, this possibility is not even parametrically satisfying, since fine tuning the 
coupling constant also implies that the vacuum energy, and so the Hubble constant during 
inflation, drops substantially. 

Turning to the radiation energy density next, our interest is to compute from the ratio 
in eq. floTS|) o.{tai)- To keep an explicit example in mind, we again take tbi = with 
C\ and C*2 as given in eqs. ( |87D and flSBl). t A i is determined by a solution once again to 
b(r) = c(r), so for a second solution to eq. ( p5|) for tai > tbi = 0. One clearly expects a 
second intersection, since b(r), which goes as (ns/3 T + 1)~ 4 ; is falling-off faster than c(r), 
which at large time goes as (k 8 / 3 t + 1)~ 2 , and after the first intersection at t B i, b(r) > c(r). 
To determine tai, only the growing mode of s a (r) is retained from (ds(r) j 'dr) 2 , to give the 
relation 



K 



8/3 



{K%/zTAI + l) 2 + K 8/3 (k 8 /3TaI + 1) - 1 = 0. 



(91) 



Solving for tai with the constraint tai > we find 



tai 



K 8/3 |_ K 8/3 



-(V2-l)-l 



so that from eq. ( pTj ) 



2(V2-1) 



-2 ' 
C 8/3 



(92) 



c{r AI ) 



K 



8/3 



16( v / 2- l) 4 



(93) 



and 



\C( T BI)J 



2^8/3- 



(94) 



Recalling from eq. (|89|) that N e e- folds of expansion requires k 8 / 3 = 1/N e , if N e ~ 50 — 70 



we find that the temperature drops by a factor 1/50 — 1/70 from the beginning to the end 
of inflation with the duration of inflation being tai — t~bi ~ N 2 ~ 2500 — 4900. 

To summarize, in this subsection we have presented an example that can be verified by 
inspection in which the scale factor expands sufficiently and then smoothly tends towards 
a radiation dominated behavior at large time. In the course of this, the temperature of the 
universe during the inflation- like stage drops between 1-2 orders of magnitude. 
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C. The Quadratic Limit 



Let us next examine the quadratic limit with the vacuum decay function 

6(r) = exp(-K 2 r). (95) 
The plots are in figure 1 for the inverse Hubble parameter 

1 _ a(r) 

H = da(r)/dr ^ 

and in figure 2 for the temperature ratio a(r) defined in eq. (|69"D. In both figures the solid 
(dashed) curve is for the vacuum decay coefficient k 2 = 0.03(0.04). The calculations were 
done by a numerical integration of the coupled scale factor equation ([R]) and stress-energy 
conservation equation (|TT|) , with the vacuum decay function in eq. (p5|) . As a cross check, 
at each iteration the so computed scale factor was substituted into the left-hand-side of 
Friedmann's equation ([8]), with the resulting energy density found from the right-hand-side 
compared with that from the numerical integration. This also means the results cross check 
separately for 6(r) and c(r) from fl66|) and (|67|) , but since this requires a(r), it is not a more 
stringent test. 

For both cases in Fig. (1), the inverse Hubble parameter starts out flat, which is charac- 
teristic of inflation-like behavior. It then veers up at a time of order r ~ 33(25) in the solid 
(dashed) case and finally tends to a slope of 2 at large times, thus becoming a radiation dom- 
inated universe on schedule. Inflation begins in the solid (dashed) case at tbi = 0.073(0.074) 
and ends at r AI = 323(228), so that the duration of inflation is tai — t bi 323(228). Dur- 
ing the inflationary period, the scale factor expands rapidly with total e-folds N e = 67(51), 
which agrees with estimates from our approximation formula eq. (|4"8"D, and the temperature 
drops by a factor (x(tai) = 11(10). 

The numerical results for tai arid ol{jai) also can be cross checked to approximate an- 
alytic expressions, Solving eq. ([85]) for the second solution at tai > Tbi, noting from 
the Appendix that K (z2(r)) dominates in the solution eq. ( |3~8"D at long time, and using 
dK (z)/dz = —Ki(z), one finds for any k 2 the general relation 

K (z 2 (r AI )) V [m) 



From |^2[ one finds that this is satisfied for z 2 ~ 1.05, so that from eq. ( j40| ) we obtain the 
approximate formula 

2 , /1.05/c 2 \ , x 

Assuming r BI <C 1 so that c{t bi ) rs 1, as for the numerical cases presented above, and using 
eq. ( ]98| ) we obtain 

a(T A i) w \V\MJT 2 . (99) 
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One can verify that the approximation formulas eqs. ( |98|) and d99|) reproduce the results for 
tai and ol(tai) respectively that were quoted above from the numerical calculations. Similar 
approximation formulas can also be obtained for the n^2 cases. 

In Fig. (2) observe the initial steep drop in a{r) for r < 1. There is a very short 
initial transient period in which the initial radiation energy density stabilizes, followed by 
a steady state stage. For both cases in Fig. (2), we started with a radiation energy density 
c(0) = 1.5. Thus initially the first term on the right-hand-side of equation ( |TTD (the "sink 
term") rapidly depletes p r (t) (equivalently c(r) in the rescaled theory) until an approximate 
balance is reached by the second term (the "source term"), after which steady state is 
reached. The initial conditions on the radiation energy density have a mild effect on the 
longtime behavior. For example, increasing c(0) by a factor of 500 has less than a one 
percent effect on N e . Without the source term, which arises from vacuum energy depletion, 
all the radiation energy would rapidly red-shift away, as in supercooled scenarios. 



V. CONCLUSION 

There are two possibly concerning or possibly predictive outcomes of scenarios occurring 
entirely in the big-bang-like regime We believe they are general features of such scenarios, 
although we do not have proof. Firstly to attain an observationally consistent expansion 
factor, it does not appear possible for the post inflation temperature Tai to be the same 
order of magnitude as that just before inflation Tbi- In supercooled scenarios this is re- 
ferred to as a perfect reheating and can be achieved by adjusting the decay width, which 
controls the reheating time period, to be sufficiently large P3UISj . In big-bang-like scenar- 
ios, for observationally sufficient expansion, we generally find that Tai is at least one order 
of magnitude below Tbi- Thus for inflationary dynamics at the grand unified scale one 
expects Taj ~ (0.1 — 0.01)Mcut- In the context of grand unified theory, this implies the 
X-boson, with Mx ~ Mqut, would not participate in post-inflation baryongenesis, although 



the lighter Higgs boson still could [23|. Moreover the picture is further altered since in the 
big-bang-like class of inflation scenarios, there would be no violent discontinuities in p r (t) 
at the end of the inflation- like stage. This implies baryongenesis could commence within 
the inflation-like stage and smoothly carry-on afterwards. One can also consider long sus- 
tained big-bang-like inflation scenarios, in which the temperature drops by a few orders 
of magnitude during the inflation-like stage. For such scenarios, studies of baryongenesis 
from out-of-equilibrium decay processes at temperatures well below M GU t m &y be useful 
|2~5fl . As a final complementary note to this concern pertaining to baryongenesis, the lower 
temperature condition implies magnetic monopole suppression works effectively. 

The second point of concern for big-bang-like inflation scenarios with not too large a 
drop in the temperature during the inflation-like stage is that they generally appear to 
have not very large upper bounds on the expansion factor with e-folds N™ ax ~ 1000 but 
jymax ^ ^qq foging typical. Since observation indicates N e > 50 — 70, this is still acceptable. 
For comparison, in the solution of the Coleman- Weinberg model in new inflation, it is found 
that N e ~ 10 7 II . In general new inflation models are reported to predict very large e-folds 
N e [IB]. As one optimistic interpretation about the small e-fold constraint for scenarios 



within the big-bang-like inflation regime, this is preferred if the universe is between nearly 
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flat and open p6 . 



In this paper we have shown that for a large range of vacuum energy density decay 
trajectories, an early universe initially in a radiation dominated stage can enter an inflation- 
like stage and finally enter back into a radiation dominated stage with the radiation energy 
density suffering no sharp alterations during this motion and with a post inflationary tem- 
perature within a range consistent with observation and theory This regime differs from the 
standard inflation regime where the radiation energy density quickly vanishes at the onset of 
inflation and then is quickly regenerated at the end of inflation in a short time reheating era. 
We reemphasize that the solutions we have found are properties of the Einstein equations, 
independent of quantum field theory. We also reiterate that in the presence of nonnegligi- 
ble radiation, one need not be restricted to familiar near equilibrium quantum field theory 
methods in searching for dynamical models. This is not to preclude conventional treatments. 
In fact, despite our emphasis on the kinematic properties of the final answer and its model 
independent origin, one should note that we motivated all our results from the conventional 
dynamical picture. 

One cannot say without further investigation what the relevance of the present results 
are. It has been established by this study that sufficiently-rapid expansion behavior is more 
general than only that found in the inflation regime. As with any generalization, there 
is always a danger that it is nothing more than a mathematical novelty offering no new 
physical insight. In the present case, this does not appear to be a correct statement. Firstly, 
in light of the new list of options, there seems no special reasons that favor the supercooled 
limit to any of the other possibilities demonstrated here. Also if the naturality principle 
carries the interpretation that any possibility not otherwise ruled out by observation nor 
theoretical common sense is a candidate solution, then again the present generalization 
has substance. Finally in conjunction with warm inflation 0, a suggestive solution to the 
amplitude fluctuation problem presents itself. However a dynamical explanation for large 
dissipation, which is needed for that scenario, requires investigation. On the other hand, 
cosmic string formation [2^-29,^|, which is typically considered a post-inflation mechanism 



for large scale structure or beginning at the end stage of a supercooled scenario |KJ, could 
be a possible mechanism within a large period of a big-bang-like inflation scenario. 

The most interesting result from this study is the finding of an inflation-like regime of 
scale factor behavior that asymptotes to the radiation dominated regime without a reheating 
stage. However returning to the introductory comments, the solutions presented here have 
applicability also to the class of supercooled inflation scenarios. In addition by decreasing 
any of the coefficients B in eqs. (|33|) and fl34|), one can smoothly interpolate from the 
big-bang-like stage of inflation to the supercooled stage of exponential expansion. There 
are many possibilities suggested by our results from mild but long sustained accelerated 
expansion to the standard exponential inflation. With any of these, due to the presence of 
radiation, the dynamic explanation may require the range from familiar methods of finite 
temperature quantum field theory to a full nonequilibrium statistical mechanical treatment 
R. Both further theoretical research and experimental information is needed to narrow the 



5 Some considerations for reheating, such as in [36], may be useful also here. In addition further 
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possibilities. 

There are two extreme points of view that one might adopt. One is that the very early 
universe mostly wants to be radiation dominated but just sneaks into a inflation-like phase 
for a little while. The other is that the very early universe has trouble containing radiation 
energy just after the initial singularity, so copiously inflates until some heating or reheating 
mechanism finally stabilizes the radiation energy. The former is the mildest modification 
of pre-inflation era thinking and the latter reflects present thinking. For now, experiment 
and theory do not indicate a strong preference for either viewpoint. However it is a useful 
exercise to view the problem from both extremes, since from either end the other looks like 
a remote limiting case. This in our opinion is symptomatic of a misunderstanding about the 
radiation energy content in the very early universe. As such, we believe theories that make 
no presumptions about the radiation content all during the early universe better represent 
the present status of experimental information about this time period. Thus, allowing for 
any of the possible scale factor behaviors derived here appears a more realistic starting point 
to further study. 
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APPENDIX A: 

This appendix contains some properties of the Modified Bessel functions K v (z) and I u (z), 
that are useful for the results in the text. At small z, the asymptotic behavior for v — are 

J„(|z|->0)~1 (Al) 
K (\z\ -> 0) ~ -m|z| (A2) 

and for v ^ 



examination could be made of the possible nonequilibrium potentials (or free energy functionals) 



that can form. An example from spinodal decomposition is [37|, which starts from a master 
equation and attempts to deduce the free energy functional. This treatment was for a conserved 
order parameter, whereas for inflaton dynamics a similar treatment is needed for a nonconserved 



order parameter. Other approaches which may be useful are in [38|. Finally the works in [39] offer 



guidance in formulating the nonequilibrium problem in an expanding universe. 
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# m (M-o)~^)(~m) 

where eq. (|A3|) is valid for all v except v ^ — 1, — 2, • ■ -. 
At large z for all z/ 

- oo) - GXP(N) 



K v (\z\ -> oo) ~ J_jLexp(-|*|). 

Finally recall that 

K u (z) = K^(z), 

and for a negative argument 

4(-|^|) = e^(|^|) 
/^(-H) = e-^(H)-i7r/ I/ (H). 
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FIGURE CAPTIONS 

figure 1: The inverse Hubble parameter for a quadratic slow-roll potential with vacuum 
decay coefficients for the solid (dashed) cases n 2 = 0.03(0.04). The initial conditions are 
c(0) = 1.5. 

figure 2: The ratio a(r) of the universe's temperature at cosmic time r to that at the 
beginning of inflation tbi for the same cases as in Fig. 1. For the solid (dashed) curve, the 
inflation-like stage begins at tbi = 0.073(0.074) and ends at tai = 323(228). In both cases 
the temperature of the universe drops by about a factor 10, with e-folds 67 (51). 
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